Bounds for weight distribution of weakly self-dual codes 

Vwani P. Roychowdhury Wid Farrokh Vatan* 
February 8, 2008 



Abstract 

Upper bounds are given for the weight distribution of binary weakly self-dual codes. To get these 
new bounds, we introduce a novel method of utilizing unitary operations on Hilbert spaces. This 
method is motivated by recent progress on quantum computing. This new approach leads to much sim- 
pler proofs for such genre of bounds on the weight distributions of certain classes of codes. Moreover, 
in some cases, our bounds are improvements on the earlier bounds. These improvements are achieved, 
either by extending the range of the weights over which the bounds apply, or by extending the class of 
codes subjected to these bounds. 

Keywords: Weight distribution, Self-dual code, Hilbert space, Unitary operation. 

1 Background 

For a random linear [n, k] code C with weight distribution (Aq, A±, . . . , A n ) it is known that the expected 
value of the normalized weight distribution, i.e., ~^A W , is the same as the normalized binomial distribution 
3tt(") (see, e.g., [14, p. 287]). So for such a code the expected value for the number of codewords of 
weight w is 2(n 1 _ fc) ( ™) . The problem to determine which explicit classes of codes have binomial weight 
distribution has been investigated in several papers (see references). Most of these results are about the 
BCH codes, their extensions, or their dual codes. For example, it is shown that for these codes the number 
of codewords of weight w is 

±_= 2^ (Q + ^) ' 
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for w in some range, and the error term E w tends to zero when n tends to infinity. 

There are also bounds for other classes of liner codes. Let C be a doubly-even self-dual [n, n/2, d] 
code C with weight distribution (Aq, Ai, . . . , A n ). Let 5 = -, then in [11] it is shown that 



A w < 2(^(f)-i)", 



(1) 



if f G [c, 1 - c], where 



1 _ 66- 1 + Vl- 85 + 325 2 

2 ~ V 8(1 - 5) ' 



(2) 



This shows that for this class of codes, the weight distribution around | is upper-bounded by the binomial 
distribution. 

In [12] an upper bound for weight distribution of the dual of extended BCH codes is given. Let 6 of 
length n = 2 m be the dual of the extended code of a t-error correcting BCH code. Then dim(C) < mt; 
i.e., |C| < n l . If (Aq, Ai, . . . , A n ) is the weight distribution of C, then for w > 



t -|j + 2 we have 



4v / 27rrari t 



1 < 

W — \2y/n(t+l)-n+2w\ 



(n — 2w) , / 1 \ \ 

" • (i + o(t))- 



(3) 



In [18] several bounds for the weight distribution of subfield subcodes of algebraic-geometric codes 
is given. This class of codes contains important classes of codes, such as the binary BCH and Goppa 
codes. For an [n, k] binary code of this type, with weight distribution (Aq, A\, . . . ,A n ), they derived the 
following bound 



1 fn 

2 n ~ k \w 



< c\n 2 



(4) 



and for the special case of w = they get the following bound 



fn— k 



n 

2 



(5) 



Here c\ and C2 are two constants, both much larger than yfe. 

In this paper, we apply a novel approach based on unitary operations on Hilbert spaces, and derive 
bounds for the weight distribution of another class of linear codes. In particular, we study the class of 
weakly self-dual codes, i.e., the class of codes S such that 6 C G 1 . We show that, for < w < ^, the 
number A w of codewords of weight w in C satisfies the following bounds 



A w < 25^(f )», 



(6) 
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and 



A w < y/e(n-w + 1)^. (7) 

If we compare our bounds with previously known bounds (1), (3), and (4), we realize that these new 
bounds, for some values of w and in the intersection of their corresponding classes of codes, give a better 
estimate than the old bounds. For example, (6) holds for any value of w, < w < and it applies also to 
the special case of doubly-even self-dual codes, and, hence, comparable to the bounds in [1 1]. The bound 
in (1) applies to doubly-even self-dual codes, but holds only in the interval [c, 1 — c], with c defined by 
(2). So if ^ < H 2 ~ 1 {\) = 0.1100 • • • then c > 0.27. Hence, for this choice of ^, (6) is a better bound 
than (1) for — < 0.27, since (1) does not even hold for these values of w; we should note, however, that 
(1) is a better bound if ^ > 0.27. One can also verify that (7) gives a better bound than (4), as the constant 
ci is much larger than yfe. 



2 Unitary operations on Hilbert spaces 

We derive our bounds via actions of unitary operations on Hilbert spaces. Toward this end, we find it 
more comfortable to use the language of "bra-kef" of quantum mechanics (see, e.g., [3]) as it is used in 
the theory of quantum computation (see, e.g., [15] and [1]). We briefly describe the necessary notions and 
notations. 

Consider the two-dimensional Hilbert space "K = C 2 . We denote its standard basis by {|0) , 1 1)}; i.e., 
|0) = (1, 0) and |1) = (0, 1). We consider also the tensor product 

<K® n = C 2 <gi C 2 <g> • • • <g> C 2 

of n copies of C 2 . So "K® n is a 2 ra -dimensional Hilbert space isomorphic with C 2 ". We represent the 
standard basis of 'H® n by 2 n products of the form 

|ci) (g> |c 2 ) ® • • • <8> |c n ) , 

where q G {0, 1}. For simplicity, we write \c\C2 ■ ■ ■ c n ) instead of \c\) <g) fa) <8> ■ ■ ■ <8> |c n ). The Euclidean 
length on l K® n is defined in the natural way; we denote the length of the vector \a) G "H® n by || \ a) \\. For 
example, 3i® 2 = C 2 (g> C 2 is a 4-dimensional Hilbert space with 

{|00),|01),|10),|11)} 

as its standard basis. If 

|o) = o|00) +6|01) + c|10) + d|ll) G :h® 2 
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then || | a) || = (aa* + bb* + cc* + dd*) 1 , where * stands for the complex conjugate. 

An m x m matrix M is unitary if Aft ■ M = I m , where I m is the identity matrix and Aft is the 
adjoint matrix of M; i.e., Aft = (M tr )*, where "tr" denotes the transpose. A linear operation on the 
m-dimensional Hilbert space C m is called a unitary operation if it is represented by a unitary matrix. 
Note that every unitary operation is a length-preserving operation. 

We denote the group ofmxm unitary matrices by U(m). The general form of a matrix in U(2) is 

/ e ^ Q+ T)cos0 e^+^sinfl^ 
\- e -i{P-i) S i n e -i(a-7) cos 6>J ' 

Suppose that f7i, G U(2), then the tensor product C/i <8> C/2 in C 2 C 2 is defined in the natural way. 
For example, if C/i(|0)) =a|0) + 6|l) andC/ 2 (|l)) = c|0) + d|l) then 

C/i®C/ 2 (|01)) = ?7i(|O))0f/ 2 (|l)) 

= (o|0) + 6|l))®(c|0) + d|l)) 

= ac |00) +ad|01) +6c|01) . 

3 Bounds for weights 

For any real number 8, consider the unitary operation Re G U(2) defined by the following matrix 



Re = 



sin cos V 
, cos 9 — sin ( 



Then the action of Rq on a vector \c), c G {0, 1}, can be written as follows: 

R e (\c})= {-^ ac (.sme) 1 ' c - a+2ac (co S e) c+a ' 2ac \a}. 

ae{0,l} 

Now let Sg = Rg® n , i.e., the tensor product of n copies of Rg. Thus Sg G U(2 n ). The following lemma 
provides a closed form for the action of Sg. 

Lemma 3.1 For any vector |c) in the standard basis, i.e., c £ {0, l} n , we have 

Sg(\c})= (-l) c ' a (sin^)"- w ' (c+a) (cos^) w '( c+a ) \a), (8) 

ae{0,l} n 

where c • a is the inner product of c and a as real vectors, wt(x) denotes the Hamming weight of a binary 
vector x, and the addition of binary vectors is considered over GF(2). 
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Proof. We first show that 

S e (\c))= ^ (-I) c - a (sin0) n - w *( c )- wt ( a ) +2c - a (cos0) w '( c ) +w *( a )- 2c - a |a), (9) 

ae{0,l}" 

We prove the identity (9) only for n = 2; the proof in the general case is quite similar. Let c = (c±, c 2 ) 
and a = (cti, G12). 

S e (\c)) = R e (\ Cl ))(E)R e (\c 2 )) 

= ( Yl (-l) aiCl (sm^) 1_Cl " ai+2aiCl (cos0) Cl+ai - 2aiCl |ai>) 

ai£{0,l} 

( 2 (-^^^(sin^) 1 -^-"^ 2 " 2 ^^^^^ 2 - 202 ^ |a 2 )) 
a2e{o,i} 

_ ^_^^aici+a2C2^ g j n ^^2-(ci+C2)-(oi+02)+2(aici+a2C2) 
ai,a2S{0,l} 

( COS 60 (ci+C2)+(ai+a2) ~ 2(aiCl+a2C2) |aia 2 ) , 

which is the of form of (9). We note that 

wt(c) + wi(a) — 2c • a = wt(c + a). 

Therefore, (9) can be rewritten as (8). I 

Now we consider a linear code C C {0, l} n of dimension k, and the corresponding unit vector 

vz cee 

By applying the unitary operation Sq on the unit vector |S) we get 

Se(\e}) = ~^=Y1 £ (-I) c - a (sin0) n - w '( c+a )(cos0) w '( c+a )|a). 
v 2k ceeae{o,i} n 

This can be rewritten as follows 

s e (\e)) = -^=Yl E( si ^) n_wt(c+a) ( co ^) wt(c+a) ) l a ) + l remamder > • 

^ 2k aee^ Vcee / 
Since ||^(|S))|| = 1, the following lemma follows. 

Lemma 3.2 For any k-dimensional linear code C of length n, and any real number 9 we have 

Yk £ (E^r^-wr^) - L (10) 

aee x Vcee / 
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Lemma 3.3 For any k-dimensional weakly self-dual linear code C of length n, i.e., C C Q- 1 , and any real 
number 6 we have 



(sin^) n - w *( c )(cos^) w * (c) 



< 2(n-2fe)/2 



(11) 



Proof. We apply Lemma III.2 to the code Q 1 . The result is 



wt(c+a) I 



Since CCC 1 , for every a G C we have 



(sin0) n - w *( c+a )(cos0) w *( c+a ) = Y (sm^) n - wt(c) (cos^) wt(c) , 



and the lemma follows. I 



Lemma 3.4 Let C be a weakly self-dual code with weight distribution (Aq, A±, . . . , A n ). Then for any 
< A < 1 we have 



n/2 

^A 2j A^ < (l + A)"/ 2 , 0<A<1. 

3=0 



(12) 



Proof. We first apply the following form of the Mac Williams identity (see [14]) to get a handier form 
of inequality (11): 

£ x n -^)y^) = _L J2(x + y) n ~ w ' (,x) (x - y) w '^. 
This way, from (11), we obtain the following inequality: 



^(sin 6 + cos 0)™- w *( c ) ( S in 6 - cos #) wt ( c ) 
cee 



< 2 n/2 . 



(13) 



Now suppose that \ < 8 < |, = sin 6* + cos0, and y 7 ?; = sin# — cos^. Thus 1 < u < 2 and 
v = 2 — u. Then inequality (13) can be written as 



(14) 



cee 
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Since S is weakly self-dual then Aj = 0, for every odd index j. Therefore, (14) can be written as follows 
(assume n is even): 

n/2 

J2 A 2j u%- j {2 - u) j < 2"/ 2 , 1< u < 2. (15) 

3=0 

Let ^=M = A, then < A < 1 and (15) became 

n/2 

^2 A ^ xj ^ (1 + A) n/2 .I 

3=0 

Theorem 3.5 For every weakly self-dual code 6 with weight distribution (Aq, A±, . . . ,A n ) we have 

n 



and 



A w < 22^(f )n < w < -. (16) 



A w < Ve{n-w+ l) w/2 , 0<w<^. (17) 



2 

Proof. For < w < § , let a = % and A w = 2^ n . From (12) it follows 

2^Af <(1 + A)§. 

Therefore, 



(3< -I a log 2 A + ilog 2 (l + A), 0<A<1. (18) 



For fixed < a < i, let 



F(A) = -ialog 2 A + ilog 2 (l + A). 



A simple calculation shows that, on the interval < A < 1, the minimum of F(X) is achieved for A = 
Thus, for < a < \, 

P<^(-a log 2 a - (1 - a) log 2 (l - a)) = ^H 2 (a). 

Therefore, we have proved (16). 

On the other hand, from (12) we have 

/ \ w/2 

Let A = ^h^, and note that (l + jf < e. This proves (17). I 
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